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1 Introduction 



In recent years, the problem of stability and stabilizability of switched systems has attracted increasing atten- 
tions (see for instance [1, 4, 5, 7, 10, 12, 13, 16]), and still many questions remain unsolved. 

In this paper, we address the problem of existence of common polynomial Lyapunov functions for linear 
switched systems. 

By a switched system, we mean a family of continuous-time dynamical systems and a rule that determines 
at each time which dynamical system is responsible of the time evolution. More precisely, let {/„ : u G U} 
(where U is a subset of E m , m G N) be a finite or infinite set of sufficiently regular vector fields on a manifold 
M, and consider the family of dynamical systems: 

x = f u (x), x e M. (1) 

The rule is given by assigning the so-called switching function, i.e. a measurable function u(.) : [0, oo[— > U C 
K m . Here, we consider the situation in which the switching function is not known a priori and represents some 
phenomenon (e.g. a disturbance) that is not possible to control. Therefore, the dynamics defined in (1) also 
fits into the framework of uncertain systems (cf. for instance [9]). 

In the sequel, we use the notations u G U to label a fixed individual system and u(.) to indicate the switching 
function. 

These kind of systems are sometimes called "n-modal systems" , "dynamical polysystems" , "polysystems" , 
"input systems" . The term "switched system" is often reserved to situations in which the set U is finite. For 
the purpose of this paper, we only require U to be measurable. For a discussion of various issues related to 
switched systems, we refer the reader to [5, 12, 13]. 

A typical problem for switched systems goes as follows. Assume that, for every u G U, the dynamical 
system x = f u (x) satisfies a given property (P). Then one can investigate conditions under which property (P) 
still holds for x = f u (t){ x )i where u(.) is an arbitrary switching function. 

In [1, 7, 10, 11], the case of switched linear systems was considered: 

x(t) = A u{t) x(t), xet", A„er x ", (2) 

where n is a positive integer and w(.) : [0, oo[— > U is a (measurable) switching function. For these systems, the 
problem of asymptotic stability of the origin, uniformly with respect to switching functions, was investigated. 

Next, we set A := {A u : u G U } and, to simplify the notation, we still call switching function the 
measurable matrix- valued map A(.) := A u i\. In this way, the switching system (2) reads: 

x(t) = A(t)x(t), where x G R™, and A(.) : [0, oo[— > A is a measurable map. (3) 

In the following, we assume that: 

(HO) the set A is a compact subset of the set ofnxn real matrices. 

Moreover, the set of switching functions, denoted by A., is the set of measurable functions A(.) : [0, oo[— > A. 
With our assumptions, for every switching function A(.) and initial condition xo G R", the corresponding 
(Caratheodory) solution of (3) is defined for every t > 0. We use 4>t ( x o) to denote the flow of (3) at time 
t > corresponding to the switching function A(.) and starting from xq. 
Let us recall usual notions of stability used for the system (3). 

Definition 1 Consider the switching system (3). We say that the origin is: 

(S) stable , if for every A(.) G Ji and e > 0, there exists 5 > such that \\4>f~^'\xo)\\ < e for every t > 0, 

Nil 

(US) uniformly stable, if it is stable with 5 not depending on A(.). 
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(U) unstable if it is not stable (i.e. if there exists A(.) £ Ji s.t. the system x(t) = A(t)x(t) is unstable as a 
linear time-varying system.) 

(AS) asymptotically stable, if it is stable and attractive (i.e. there exists 8' > so that, for every A(.) £ Ji 
and x £ R" with \\xq\\ < 6', we have lim^oo \\<f>f (xq)|| = 0). 

(UAS) uniformly asymptotically stable if it is uniformly stable and if, for every e' > and 8' > 0, there exists 
T > such that for every switching function A(.) £ Ji, t >T and \\xq\\ < 8' , we have \\4>f (^o)!! < 

(GUES) globally uniformly exponentially stable, if there exist positive constants M, A such that: \\(f>f^\xQ)\\ < 
M e- At |jx || ; for every x £ R™, t > 0, A{.) £ A. 

Due to the fact that the dynamics is linear in the state variable, the local and global notions of stability are 
equivalent. More precisely, it was proved in [3] that, for system (3) subject to HO, the three notions AS, UAS, 
GUES and the notion of attractivity are all equivalent (see also [10, 13]). In addition, if the system is unstable, 
then there exists a switching function ^4(.) £ Ji and an initial condition xq such that lim^oo ||<^^'^(a;o)|| — ► oo. 
In the following, we just refer to the notions of stability, instability and GUES. 

Remark 1 Since for the stability issue, a system of type (3), subject to HO, is uniquely determined by a 
compact set A of n x n real matrices, we identify A with the corresponding system for the rest of the paper. 
For instance, when we say that A is GUES, we mean that the corresponding system of type (3) is GUES. 
We will often consider the problem of determining whether a system, belonging to a certain class C of systems 
of type (3) subject to HO, is GUES or not. Notice that fixing such a class of systems means to fix a set of 
compact subsets of R nx ™ i.e. C can be identified with a subset of {A C R nxn : A compact}. 

For a system (3) subject to HO, it is well known that the GUES property is a consequence of the existence 
of a common Lyapunov function. 

Definition 2 A common Lyapunov function ( CLF for short) V : 1™ — > R + ; for a switched system (S) of 
the type (3), is a continuous function such V is positive definite (i.e. V(x) > 0, Va; ^ 0, V(0) — 0) and V is 
strictly decreasing along nonconstant trajectories of(S). 

Vice- versa, it is proved in [10] that, given a GUES system of the type (3) subject to (HO), it is always possible 
to build a C°° common Lyapunov function. 

Anyway, the problem of finding a CLF or proving the nonexistence of a CLF is in general a difficult task. 
Sometimes, it is even easier to prove directly that a system is GUES or unstable. An example is provided 
below by bidimensional switched systems. 

1.1 Single-Input Bidimensional Switched Systems 

Consider a single input bidimensional system of the type: 

x(t) = u(t)Ax(t) + (1 - u(t))Bx(t), (4) 

where a; £ R 2 , A and B are two 2x2 real Hurwitz matrices and u(.) is a measurable function defined on 
R + and taking values in U equal either to [0, 1] or {0, 1}. In the sequel, we call S the class of bidimensional 
systems of the above form. This class is parameterized by couples of 2 x 2 real Hurwitz matrices. 

Remark 2 Whether systems of type (4) are GUES or not is independent on the specific choice U — [0, 1] or 
U = {0, 1}. In fact, this is a particular instance of a more general result stating that the stability properties of 
systems (3) subject to HO only depend on the convex hull of the set A, see Proposition 1 below and Appendix 
B. 
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In [15], the authors provide a necessary and sufficient condition on the pair (A, B) to share a quadratic 
CLF, but Dayawansa and Martin showed in [10] that there exist GUES linear bidimensional systems not 
admitting quadratic CLF. They also posed the problem of finding the minimal degree of a polynomial CLF. 
More precisely, the problem posed by Dayawansa and Martin is the following: 

Problem P: Define 'E.gues C 3 as the set of GUES systems of the type (4). Find the minimal integer 
m such that every system of Egues admits a polynomial CLF of degree less or equal than m. 

Remark 3 In the problem posed by Dayawansa and Martin, it is implicitly assumed that a GUES system 
always admits a polynomial Lyapunov function and one of our results (our Theorem 1 below) indeed confirms 
that fact. 

As for the GUES issue, it was completely resolved in [7], where a necessary and sufficient condition for a 
system of type (4) to be GUES was found directly, without looking for a CLF (see Section 3 and Appendix 
A for more details). This is a typical example in which it is easier to study directly the stability rather than 
looking for a CLF. 

1.2 Sets of Functions Sufficient to Check GUES 

The concept of Lyapunov function is useful for practical purposes when one can prove that, for a certain class 
of systems, if a CLF exists, then it is possible to find one of a certain type and possibly as simple as possible 
(e.g. polynomial with a bound on the degree, piecewise quadratic etc.). 

More precisely, consider a class C of systems of type (3) in I™ subject to HO, in the sense of Remark 1. 
One would like to find a class of functions Sc, identified by a finite number of parameters, which is sufficient to 
check GUES for systems belonging to C i.e., if a system of C admits a CLF, then it admits one in Sc- Once 
such a class of functions is identified, then in order to verify GUES one could use numerical algorithms to 
check (by varying the parameters) whether a CLF exists (in which case the system is GUES) or not (meaning 
that the system is not GUES). 

For instance, a remarkable result for a given class C of systems in R n could be the following: 

Claim: there exists a positive integer m (depending on n) such that, whenever a system of C admits a 
CLF, then it admits one that is polynomial of degree less than or equal to m. In other words, the class of 
polynomials of degree at most m is sufficient to check GUES for the class C. 

If this result were true, one could use numerical algorithm to check, among all polynomial of degree m (varying 
the coefficients), if there is one that is a CLF. Unfortunately, this claim is not true , even for the simplest non 
trivial case of class of systems in R 2 , namely systems of type 5 (cf. Equation (4)). 

The next definition formalizes the idea of class of functions sufficient to check GUES. 

Definition 3 We say that a subset S o/C°(R n ,R) is finitely (or m- finitely) parameterized if there exist il C R m 
for a positive integer m and a bijective map $:!1C R m — > S c C°(R™,R). A subset S o/C°(R™,R) is said 
to be sufficient to check GUES for a class C of systems of type (3) in R" (SSF for short), if every GUES 
system of C admits a CLF in S. A subsets ofC (R n ,R) is said to be an (m-parameters) finite set of functions 
sufficient to check GUES for a class C of systems of type (3) in W 1 (finite-SSF for short), if S is m-finitely 
parameterized and is an SSF for C . 

If a subset S of C °(R n ,R) is not finitely parameterizable but is an SSF for a class of systems C , we call S 
an oo-SSF. 

Remark 4 In [6] , a concept similar to those introduced in the previous definition was provided and it was called 
"universal class of Lyapunov functions" . 

Using the previous definitions, the results and the problem formulated in [10] can be rephrased in the following 
way: 

Rl for systems (3) subject to HO, the set C°°(R",R) is an oc-SSF; 
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R2 for linear bidimcnsional systems of the class 3 (cf. Formula (4)), the set of quadratic functions is not a 
SSF; 

P let V m be the set of polynomial functions of two variables with degree at most to. What is the minimal m 
such that V m is a finitc-SSF for the linear bidimcnsional systems of the class S? 

Remark 5 Notice that, to check numerically the existence of a CLF using the concept of finitc-SSF, one needs 
some regularity properties of the functions of the family, with respect to the parameters (at least continuity). 
Anyway, this discussion is out of the purpose of this paper. 

1.3 Main Results 

We first prove that the implicit assumption of Dayawansa and Martin (i.e. that a linear GUES switched 
system always admits a polynomial CLF, cf. Remark 3) is correct in R" (and in particular for bidimcnsional 
systems of type (4)). 

Theorem 1 If the origin is a GUES equilibrium for the switched system (3) subject to HO, then there exists 
a polynomial CLF. 

The above result can be stated equivalently as follows. 

Theorem Ibis The set of polynomials from R™ to R, is an oo-SSF for linear switched systems (3) sub- 
ject to HO. 

The proof of Theorem 1 is given in Section 2 and the starting point is the construction of a homogeneous 
and convex CLF W, following the corresponding argument of [10]. The main idea is then to seek for a (homo- 
geneous) polynomial W whose level sets approximate, in some suitable sense, those of W and, finally to show 
that W is also a CLF. 

A related result has also been obtained in [6] , using an intermediate approximation with polyhedral functions 
(this step is implicit in our proof) and starting from the case of discrete approximating systems (so-called Eulcr 
approximating systems). 

The core of the paper consists of showing that problem P does not have a solution, i.e. the minimum de- 
gree of a polynomial CLF cannot be uniformly bounded over the set of all GUES systems of the form (4). 
More precisely, we have the following: 

Theorem 2 Let 'E.gues C H be the set of all GUES systems of the type (4). If (A, B) is a pair of 2 x 2 real 
matrices giving rise to a system ofBcuES, let m(A, B) be the minimum value of the degree of any polynomial 
CLF associated to that system. Then m(A,B) cannot be bounded uniformly over "E-gues ■ 

The above result can be stated equivalently as follows. 

Theorem 2bis Let V m the set of polynomial functions from R 2 — > R of degree at most to. Then V m is 
not a finite-SSF for S. 

The proof, given in Section 4, is based on ideas developed in [7], where necessary and sufficient conditions 
for GUES of systems (4) are provided. We build a sequence of GUES systems corresponding to a sequence 
of pairs of matrices (Ai,Bi), i > 1. The sequence of systems is chosen in such a way that the limit system 
is uniformly stable but not attractive. In particular, that limit system admits a nontrivial periodic trajectory 
whose support T is a C 1 but not a C 2 submanifold of the plane. To each GUES system of the sequence, one 
considers any polynomial CLF g whose degree is at most to. We prove that a subsequence of (V^. jj.) 
converges to a non zero polynomial function V (of degree at most to) which admits T as a level set. Since T is 
not analytic, a contradiction is reached. 

Remark 6 The result given by Theorem 2 generalizes to dimensions higher than 2 as follows. Let (Ai, Bi), i> 
1, be a sequence of 2 x 2 matrices such that i) the corresponding systems of type (4) arc GUES, ii) the limit 
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is uniformly stable but not attractive. As explained above, for this sequence of systems it is not possible to 
build a sequence of polynomial CLF of uniformly bounded degree. Consider now the sequence of systems in 
R™, n > 2, of the form x = uAix + (1 — u)BiX corresponding to the matrices: 
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Each system of the sequence is GUES but the limit system is not (it is just uniformly stable). Now, if V^. 
i > 1, are the corresponding polynomial CLFs, then they cannot be polynomials of uniformly bounded degree 
since this is not true for the restriction of V^. g. to the first two variables. 

Remark 7 (extension to piecewise polynomial functions (PPF)). Another class of functions commonly 
used to check GUES is that of piecewise quadratic functions or more generally piecewise polynomial functions 
(PPF for short). Here, by a PPF, we mean a continuous function V € C°(R", R) together with a finite number 
q of cones Kj, 1 < j < q, based at zero and partitioning R™ so that V is a polynomial function of degree dj on 
Kj, 1 < j < q. We refer to m(V) := max{g,<ii, ...,d q } as the total degree of V. 

It is tempting to state a version of problem P by replacing polynomial functions of degree at most m with 
PPFs of total degree at most m. 

Again, the PPF version of problem P does not have a solution for n = 2, i.e. the minimum total degree of 
a piecewise polynomial CLF cannot be uniformly bounded over the set of all GUES system of the form (4). 
The argument is a simple extension of the proof of Theorem 2 and it is briefly mentioned in Remark 13. 

The last results of the paper concern the existence and the characterization of a finite-SSF for systems of 
the type (3) subject to HO. 

Let us define the convex scmicone generated by a set D C R ra as the set of points Xx with A > and 
x G co(D), where co{D) denotes the convex hull of the set D. With this definition, the point x = does not 
belong to the convex semicone generated by a set D, if ^ co(D). 

First of all, we prove the following (see Appendix B for the argument). 

Proposition 1 For every compact subset A ofW lxn (i.e. verifying HOj, let Sj^ be the system of the type (3) 
associated to A. Then, for A and A verifying HO and generating the same convex semicone, Sj^ is GUES 
(resp. uniformly stable) if and only ifSj^' is GUES (resp. uniformly stable). 

Based on converse Lyapunov theorems, one can deduce some trivial existence results for finite-SSFs. For 
instance, consider a class C of systems of type (3) in R™ subject to HO and satisfying the following property: 
for every A E C, the convex hull of A is generated by at most k matrices n x n, where k is a positive integer. 
One can build a finitc-SSF for the class C as follows. If C does not admit any GUES system, the positive 
definite function x i— > ||a;|| 2 will detect the instability of every element of C. In this case, the subset of C°(R™, R) 
reduced to x ||x|| 2 is a finite-SSF for C. Otherwise, C admits at least one GUES system. Assume first 
that a finite-SSF S can be provided for the GUES systems of C and let V & S be a CLF for a fixed GUES 
system of C. By associating V to every unstable system of C, it is clear that S becomes a finite-SSF for the 
whole class C. 

Therefore, we may simply assume that C is made of GUES systems. Thanks to Proposition 1, the class 
C can be parameterized by fc-tuples of n x n matrices, defined up to their norm. In this way, a k(n 2 — 1))- 
parameters finite-SSF is provided for the class C . For instance, the class 5 of two-dimensional systems of type 
(4) admits a 6-parameters finite-SSF. 

The above construction is not explicit and therefore is not useful to check GUES. Similarly to Lyapunov 
functions, it is then clear that the real challenge for finite-SSFs concerns their explicit characterization. For 
classes of systems of type (3) in R™, n > 3, that issue is completely open in general. In dimension two, using 
the necessary and sufficient conditions for GUES given in [7], we provide an explicit 5-parametcrs finite-SSF 
for S. This is the content of Section 5. 
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Clearly, S can be parameterized by the pairs (A, B) of 2 x 2 real Hurwitz matrices, where both A and B 
are defined up to their norm. The construction of the explicit finite-SSF goes as follows. As done previously, 
we may assume that the pair (A, B) gives rise to a GUES system. By taking advantage of the complete 
characterization of GUES systems of the class 5 given in [7], one can explicitly associate to every GUES pair 
(A, B) a CLF as explained next. We start by defining, from (A, B), a pair (A, B) giving rise to a system of 5 
which is uniformly stable but not attractive. Such a system admits a closed trajectory whose support T is a 
simple Jordan closed curve (cf. Sections 3, 4). We then construct a homogeneous positive definite function V 
whose level set 1 is T. We finally show that V is a CLF for (A, B). Since the set of (A, B) built from the GUES 
pairs {A, B) can be parameterized by using five parameters, we end up with a five-parameters finite-SSF. 

1.4 Structure of the Paper 

Section 2 contains the proof of Theorem 1. In Section 3, we recall the main ideas from [7] needed for the rest 
of the paper. For sake of completeness, we provide in Appendix A the full statement of the main result of [7]. 
The proof of Theorem 2 is given in Section 4 and the explicit construction of a five-parameters finite-SSF for 
systems of type (4) is provided in Section 5. Finally, in Appendix B, we prove Proposition 1. 

2 Existence of Common Polynomial Lyapunov Functions 

In this section, we prove Theorem 1. The starting point of the argument follows the first part of the proof of 
an analogous result in [10]. 

We define the function V : M" -> M+ by: 

r+oo 

V(x)= sup / ||<^«(a;)|| 2 dt . 
A(.)eAJo 

The function V is well defined since there exist positive constants C, ijl such that, for alH > and x £ K": 

Ut { \x)\\<Ce-^\\x\\. 

Note that V is homogeneous of degree 2 and continuous. In addition, we next show that V is strictly convex. 
That fact will be crucial later in the argument. Fix x, y £ K™ and x ^ y. Let A(.) be a switching function. The 
function x i— > ||0^'^(a;)|| 2 is strictly convex. Moreover, for every A e]0, 1[, by compactness of A , the expression: 

A||^ W (*)|| a + (1 - A) ||<^%)|| 2 - ||^ ( - ) (Ax+ (1 - A)y)|| 2 , 

is nonncgative for every t > and is bounded from below by a positive constant on some interval [0, t\, uniformly 
with respect to A(.). 

Therefore, dividing the integration interval into the two intervals [0, t\ and [i, +00] and taking a maximizing 
sequence of switching functions for V(Xx + (1 — A) y), we have: 

V(Xx + (l-X)y)<XV(x) + {l-X)V(y), VAe]0,l[, Vx,y€R n . 

It is shown in [10] that V is a CLF. Nevertheless, we need to consider at least C 1 Lyapunov functions, therefore 
we define: 

V(x) = [ f(R)V{Rx)dR. 

JSO(n) 

where / : SO(n) — > [0, +00 [ is a smooth function with support on a small neighborhood of the identity matrix 
and/ so( „ ) /( J R)d J R=l . 

In [10], it is also shown that V is a smooth CLF except at the origin. Moreover, since V is homogeneous of 
degree 2 and strictly convex, it follows that V also satisfies such properties. 
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We consider now the function W(x) = yV(x) , which is a continuous, positively homogeneous CLF. Therefore, 

W^ X {1) is a compact set. Using the fact that the set {x : W(x) < 1} is strictly convex, we construct a 
polynomial CLF W by approximating the level sets of W. For this purpose, we need the following preliminary 
result which describes a continuity property of the function \7W(y) ■ Dx with respect to x, y, D. 

Lemma 1 Let us set: 

M : = min [ - VW{x) ■ Dx] . 
xew-^i) DeA 

Then, for every e £ (0, M), there exists S £ (0, 1) such that, for every x, y £ W (1) with VW(y) ■ x > 1 — 6 
and every D £ A, one has: 

VW{y) ■ Dx < -e. 

Proof of the Lemma. First of all, notice that M is well defined since it is the infimum of a continuous 
function over a compact set. Moreover, M > because W is a CLF. 
Since, by homogeneity, VW(y) ■ y = W{y) = 1, we have: 

VW(y)-x = l-VW(y)-(y-x), 

and then the hypothesis is equivalent to VVF(y) ■ (y — x) < S . 

Reasoning by contradiction, assume that there exists a sequence (xj,yj, Dj) such that VW(j/j) • DjXj > —e 
and VW(yj) ■ (j/j — Xj) converges to as j goes to infinity. By compactness, we can find a subsequence of 
(xj,yj,Dj) converging to (x,y,D) and therefore, by continuity, VW r (y) • Dx > —e and \7W(y) ■ (y — x) = 0. 
Therefore y — x belongs to the tangent space at y of the strictly convex set W /_1 ([0, 1]). Since x also belongs 
to the boundary of that set, it must be y = x. It implies \7W(y) ■ Dx = WW(x) ■ Dx < —M and we reach a 
contradiction. □ 

Remark 8 Taking —x instead of x, one obtains that for every x,y £ and every DeA, then VW(y)-x < 

-l + S VW(y) ■ Dx > e. 

To conclude the proof of the theorem, we take 8 € (0, 1) corresponding to some e as in the lemma above, and 
for every y £ VF _1 (1) we consider the open sets B y = {.t <G W 1 : \7W(y) ■ x > 1 — S/2\. Since y <G B y , wc 
have that {Byjy^yy- 1 ^) is an open covering of the compact set W /_1 (l) , and therefore we can find yi, . . . , y^ 
points of 1F _1 (1) such that the union of B yk , k — 1, . . . , N, covers . 
Let us define: 

N 
k=l 

We claim that, for an integer p large enough, W is a polynomial CLF. For D £ A and x £ M. n , x ^ 0, we have: 

jv 

■ Dx = 2pJ2 {VW(y k ) ■ x) 2p -\w{y k ) ■ Dx, (6) 
fe=i 

and we want to show that \7W(x) ■ Dx < 0. By homogeneity, it is enough to do it for x £ W~ 1 (l). Set: 

K := max VW(y) ■ Dx. 

i, ! /eB'- 1 (i), DeA 

If, for some index k in {1, iV}, one has |Viy(j/fe) • x\ < 1 — S. Then: 

\{VW(y k ) ■ x) 2p -\w(y k ) ■ Dx\ < (1 - Sf^K. 

Otherwise, if the inequalities 1 — 8/2 > \S7W(yk) ■ x\ > 1 — 6 hold, then, by the previous lemma and remark, 
one has that the corresponding term in the summation must be negative. 

Finally, since by the definition of the points y k , there exist at least two distinct indices fci and &2 such that 
x £ B„, and —x £ B,,, we have that: 

{VW(y ki ) ■ x) 2p -\w(y kt ) ■ Dx < -(1 - S^^e. 
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Summing up, we deduce that: 




1 



K{N -2) (1-5 \2p-i 



) 



2e VI -8/2) 



For p large enough, the right-hand side of previous expression is negative, uniformly with respect to D G A 
and x G The theorem is proved. 

Remark 9 One can also check that the level set VF _1 (1) approximates, as p tends to +00, the corresponding 
level set of the function ma,Xk=i....,N \VW(yk) ■ x\ (which is a polytopc) and, therefore, the latter is a CLF as 



3 Necessary and Sufficient Conditions for GUES of Bidimensional 
Systems 

Consider the following property: 
(V) The bi-dimcnsional switched system given by: 



is GUES at the origin. 

In this section, together with Appendix A, we recall the main ideas from [7], to get a necessary and sufficient 
condition on A and B under which (V) holds, or under which we have at least uniform stability. The full 
statement of the Theorem is reported in Appendix A.l. 

Remark 10 Recall that, by Proposition 1 (proved in Appendix (B)), the necessary and sufficient condition for 
stability of the system (7) are the same if we assume u(.) taking values in {0, 1} or in [0, 1], or if we multiply 
A and B by two arbitrary positive constants. 

Set M{u) := uA + (1 — u)B, u G [0, 1]. In the class of constant functions the asymptotic stability of the origin 
of the system (7) occurs if and only if the matrix M(u) has eigenvalues with strictly negative real part for each 
u G [0, 1]. So this is a necessary condition for GUES. On the other hand it is know that if [A, B] = then the 
system (7) is GUES. So, in what follows, we always assume the conditions: 

HI: Let Ai,A2 (rcsp. A3,A4) be the eigenvalues of A (resp. B). Then i?e(Ai), i?e(A2), i?e(As), i?e(A4) < 0. 
H2: [A, B] 7^ (that implies that neither A nor B is proportional to the identity). 
For simplicity we will also assume: 

H3: A and B are diagonalizable in C (notice that if H2 and H3 hold then Ai 7^ A2, A3 7^ A4). 

H4: Let Vi, V 2 G CP 1 (rcsp. V 3 , V 4 G CP 1 ) be the eigenvectors of A (rcsp. B). Then V, 7^ Vj* for i G {1, 2}, 
j G {3,4} (notice that, from H2 and H3, the V^s arc uniquely defined, Vi 7^ V2 and V3 7^ V4, and H4 
can be violated only when both A and B have real eigenvalues). 

All the other cases in which HI and H2 hold arc the following: 

• A or B are not diagonalizable. This case (in which (V) can be true or false) can be treated with techniques 
entirely similar to the ones of [7]. 

• A or B are diagonalizable, but one eigenvector of A coincides with one eigenvector of B. In this case, 
using arguments similar to those of [7], it possible to conclude that (V) is true. 



well (cf. [6]). 



x(t) 



u(t)Ax(t) + (1 - u{t))Bx(t), where u(.) : [0, oo[-> [0, 1], 



(7) 
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We will call respectively (CC) the case where both matrices have non-real eigenvalues, (RR) the case where 
both matrices have real eigenvalues and finally (RC) the case where one matrix has real eigenvalues and and 
the other non-real eigenvalues. 

Theorem 3, reported in Appendix A.l, gives necessary and sufficient conditions for the stability of the 
system (7) in terms of three (coordinates invariant) parameters given below in Definition 4. The first two 
parameters, pa and pb, depend on the eigenvalues of A and B respectively, and the third parameter K. depends 
on Tr(AB), which is a Killing-type pseudo-scalar product in the space of 2 x 2 matrices. As explained in [7], 
the parameter K, contains the inter-relation between the two systems x = Ax and x = Bx, and it has a precise 
geometric meaning. It is in 1-1 correspondence with the cross ratio of the four points in the projective line 
CP 1 that corresponds to the four eigenvectors of A and B. 

Definition 4 Let A and B be two 2x2 real matrices and suppose that HI. H2. H3 and H4 hold. Moreover 
choose the labels (i) and (2) (resp. (3) and (4)) so that | A2 1 > |Ai| (resp. | A4 1 > | A3 1 ^ if they are real or 
Im{\2) < (resp. iro(A4) < 0) if they are complex. Define: 

\i + X-2 .A3 + A4 „ Tr(AB) - ±Tr(A)Tr(B) 

PA-=-i~, r-; PB ■= -1-: r-; /C := 2- 



Ai — A2 ' A3 — A4 ' (Ai — A2XA3 — A4) 

Moreover, define the following function of pa, pB,fc-' 

V:=K 2 + 2pAPBK,-{\ + p 2 A + p B ). (8) 



Notice that pa is a positive real number if and only if A has non-real eigenvalues and pa G iR, Pa/? > 1 if and 
only if A has real eigenvalues. The same holds for B. Moreover Del. 

Remark 11 Under hypotheses HI to H4, using a suitable 3-parameter changes of coordinates, it is always 
possible to put the matrices A and B, up the their norm, (cf. Remark 10) in the normal forms given in 
Appendix A. 2, where pa,Pb,K- appear explicitly (see [7] for more details). 

The parameter /C contains important information about the matrices A and B. They arc stated in the following 
Proposition that can be easily proved using the normal forms given in Appendix A. 2. 

Proposition 2 Let A and B be as in Definition 4- Then: i) if A and B have both complex eigenvalues, then 
K S 1 and \K.\ > 1; ii) if A and B have both real eigenvalues, then K, G K \ {±1}; iii) A and B have one 
complex and the other real eigenvalues if and only if K, G i$L. 

Theorem 3, stated in Appendix A.l, is the main result of [7], and gives necessary and sufficient conditions for 
(V) holding true. We next describe the main idea of the proof. All details can be found in [7]. 

We build the "worst trajectory" j Xo i.e. the trajectory (based at xo) having the following property. At each 
time t, 7x (i) forms the smallest angle (in absolute value) with the (exiting) radial direction (Figure 1 A). 

Then the system (7) is GUES if and only if, for each xo G K 2 , the "worst trajectory" j Xo tends to the 
origin. The worst trajectory is constructed as follows. We study the locus Q _1 (0) (where Q(x) :— det(Ax, Bx)) 
where the two vector fields Ax and Bx are collinear. The quantity T>, defined in Definition 4, is proportional 
to the discriminant of the quadratic form Q. We have several cases: 

• If <2 -1 (0) contains only the origin then, in the (CC) and (RC) case, one vector field points always on 
the same side of the other and the worst trajectory is a trajectory of a fixed vector field (either Ax or 
Bx). In that case, the system is GUES (case (CC.l) and (RC.l) of Theorem 3), see Figure 1, case B. 
The situation is similar in case (RR.l) (the worst trajectory tends to the origin). 

• If Q _1 (0) does not contain only the origin then it is the union of two lines passing through the origin (since 
Q is a quadratic form). If at each point of Q _1 (0), the two vector fields have opposite direction, then 
there exists a trajectory going to infinity corresponding to a constant switching function (see Figure 1, 
case C). This correspond to cases (CC.2.1), (RC.2.1) and (RR.2.1) of Theorem 3. In that situation, 
there exists u G [0, 1] such that the matrix M(u) := uA + (1 — u)B, u E [0, 1] admits an eigenvalue 
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Figure 1: Proof of the stability conditions 



with positive real part. If at each point of Q _1 (0), the two vector fields have the same direction, then 
the system is GUES if and only if the worst trajectory turns around the origin and after one turn the 
distance from the origin is decreased, (see Figure 1, cases D and E). The quantities p cci P RC , P RR 
defined in Theorem 3 (for the three cases (CC), (RC), (RR) resp.) represent the distance from the 
origin of the worst trajectory (that at time zero is at distance 1), after one half turn. This correspond to 
cases (CC.2.2), (RC.2.2) and (RR.2.2) of Theorem 3. 

• Finally (CC.3), (RC.3) and (RR.3) are the degenerate cases in which the two straight lines coincide. 

4 Non Existence of a Uniform Bound on the Minimal Degree of 
Polynomial Lyapunov Functions 

In this section, we prove Theorem 2. The starting point of the argument is to consider a pair of matrices A 
and B having both non real eigenvalues ((CC) case) and satisfying: 

V>0, K>1, p cc = l. (9) 

Such a pair exists. Indeed, Figure 2 translates graphically the contents of Theorem 3 for a fixed K, > 1, in 
the region of the (pa, Pb) -plane where pa, Pb > 0. The open shadowed region corresponds to values of the 
parameters pa,Pb for which the system is GUES. We denote by S + the open subset of the shadowed region 
where T> > 0. The curve C represents the limit case where p cc = 1. To each internal point of that curve, it is 
associated a system verifying (9), since T> > 0. A system corresponding to such a limit case is not asymptotically 
stable but just stable. Moreover, the worst trajectory is a periodic curve, whose support is of class C 1 but not 
of class C 2 (recall that the switchings occur on Q _1 (0), i.e. when the linear vector fields corresponding to A 
and B are parallel). 

Fix a point (pa, Pb) G C corresponding to (A, B). From the picture, it is clear that there exists a sequence 
of points (pAkiPBk) £ S + 1 f° r k>l, converging to (pa,Pb)- This exactly means that there exists a sequence 
of GUES pairs (Ak,Bk), k>l, such that (Ak,Bk) tends to (A, B) as k goes to oo. 

Let x = {x\,x%). For every k > 1, consider a polynomial CLF Vk = Ei<i+j<m t a ij' ] x \ x 2 °f degree at 
most mfc. Arguing by contradiction, we assume that the sequence (rrik) is bounded by a positive integer m. 
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CaseK>l 



(D>0) 



D<0 



Unstable q (Uniformly stable, but not attractive) 



Figure 2: GUES property in the space of parameters and explicit construction of a 5-parameters SSF for 
systems of type (4). S + is the region in the (pA, ps)-plane in which the system is GUES and T> > 0. 



Up to multiplication by a constant, we can choose Xii<i+j<m t \ a ij I = 1- By compactness, there exists a 
subsequence of (14) (still denoted by (14-)) which converges (uniformly on compact subsets of R 2 ) to some 
non-zero polynomial V with degree at most m. Note that V(0) = since the T4's are CLFs. 

Fix io £ I 2 , io / 0. Let T > be the period of the worst trajectory r y XQ corresponding to the pair 
(A, B), and starting at xo . Note that T is independent of xq. The curve 7^ : [0, T] — > R 2 can be seen as the 
concatenation of at most five arcs of integral curves of x = Ax and x = B x (see Figure 3) and satisfies the 
Cauchy problem: 

( x = C(t)x, 
\ x(0) = x , 

where C(t) is equal to A or B on subintervals of [0, T]. 
For k > 1, consider the Cauchy problem: 

( ± = C k (t)x, 
\ x(0) = x , 

where Ck(t) = Ak if C(t) = A and Ck(t) = Bk if C(t) = B. Then, 7^ is a trajectory of the switched system of 
the type (4) associated to (Ak, Bk). Since, the right-hand side of the previous equation is Lipschitz continuous 
in x and pieccwise continuous in t, then the solutions 7& converge uniformly to 7 Xo on [0, T]. 
We next show that V remains constant on j Xo . For k > 1 and t £ [0,T], one has: 

\\V k o lh (t) - V o 7ao (t)|| < \\V k o 7fc (t) - V o 7fc(t)|| + || V o lk (t) - V o 7 ao (t)||. 

By uniform convergence of Ufc to V and of 7 fe to j Xo , and by continuity of V, we deduce that Vfc0 7/c(t) converges 
to V o 7 Xo (t) for every fixed t. 

Since, for every k > 1, Vfe is a CLF for the switched system of the type (4) associated to (Ak, Bk), then Vfc o 7fc 
is a decreasing function and, hence, V o 7xo is non-increasing. Moreover U o 7xo (T) = 7^ o (0). Therefore, 
V oj Xg must be constant. It implies that there exists t\ > such that either V(e At Xo) or U(e s *a;o) is constant 
on [0, ti]. With no loss of generality, assume the first alternative. Since the map 1 1— ► U(e j4 *xo) is real analytic, 
it follows that V(e At xo) is constant over the whole real line. By letting t go to +00, since e At xo — > 0, we deduce 
that V(xq) = V(0) = 0. Since xq is an arbitrary non zero point of M 2 , we get that V = 0, which is not possible. 
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Figure 3: The "worst trajectory' 



Remark 12 The construction of the sequence (Ai,Bi) with unbounded degree for polynomial CLF was per- 
formed for matrices having both non real eigenvalues (that corresponds to the (CC) case). The same con- 
struction can be reproduced for the (RC) and (RR) cases. 

Remark 13 . For the PPF case (see Remark 7), the above argument can be easily modified to get that 
the minimum total degree of a piecewise polynomial CLF cannot be uniformly bounded over the set of all 
GUES systems of the form 5 (cf. Formula 4). Indeed, let T4 be the sequence of PPFs taking the value 
V£(x, y) = J2i<i+j< m a iji x i x i m tne c °ne Kf, for 1 < I < m. Here, to simplify the notation, we assume 
without loss of generality, that, for each element of the sequence, the number of cones and the degree of 
Vl(x,y) is always to. 

Each cone can be identified by a couple of angles with the x-dircction. Therefore to each function 14 we can 
associate a m-uple of angles (aj, . . . , a^,) such that the cone Kf coincides with the region between the lines 
corresponding to af and af +1 . In particular, up to subsequences, we can assume that the numbers af converge 

to ai. Similarly to the case above, we can normalize the coefficients of the CLFs VJt by YliLi Ei<i+3<m \ a ijl I = ^ 
and consider a subsequence of the coefficients converging to a^i. Then, if we define V as the PPF such that 
V(x, y) = V (x, y) = J2i<i+j< m a ijl x \ x °'i 011 the cone Ki defined by the angles ai and aj+i, it is easy to verify 
that Vk(x, y) converges uniformly on compact subsets of R 2 to V(x, y). We can conclude the proof as before 
showing that Vi(xq) = V(0) = for arbitrary xq, which leads to a contradiction. 



5 Explicit Construction of a finite-SSF for Systems of Type (4) 



In this section, we provide a 5-parameters finite-SSF for the class S of bidimensional systems of type (4). 
Recall that for what concern the stability issue, S can be parameterized by the 6-parameters family provided 
by the pairs (A, B) (of 2x2 matrices) defined up to their norm. 

As explained in the introduction, it is enough to construct a CLF for a pair (A, B) giving rise to a GUES 
system of B. We only treat the (CC.2.2) case since, in all the other cases, the construction is entirely similar. 

In the (CC) case, after a three-parameters change of coordinates, the normal form for the pair (A, B) is 
given by (cf. Appendix A. 2): 



Moreover, in the (CC.2.2) case, we have K. > 1, V > and p cc < 1, where JC := l/2(E + 1/E), V and p cc 
being respectively defined in (4) and (11). Recall that, for fixed K. > 1 (i.c fixed E > 1), Figure 2 describes, in 
the (pa, ps)-plane, the status of each point with respect to the GUES issue. 




E > 0. 
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Wc now associate, to every GUES pair (A, B), a pair (A,B) corresponding to a system of the type (4) 
uniformly stable but not attractive. Consider in Figure 2 the line segment joining the point (0,0) to (pa,Pb) 
in the S + region. That segment intersects the curve C in a point (pa,Pb)- That results from the Jordan 
separation theorem and the fact that C connects the points (VJC 2 — 1, 0) and (0, \JK? — 1). Therefore, there 
exists a ( g (0, 1) such that, for the system given by: 

A = ( ~E -if)' *-(-? -t)- ~>*=^ 

the worst trajectories r ) Xo are closed curves, i.e. P cc = 1. 
Moreover, one can easily compute: 

det(Ax, Ax) = p A (l - Q ( Xl 2 E + ^\ > , 

det(Bx,Bx) = p B (l ~ C)M 2 > . 

Therefore the vector fields Ax, Bx point inside the area delimited by a fixed worst trajectory (that is closed 
curve) of the modified switched system and so, passing to angular coordinates, the function: 

V(r,a) = ^-r, (10) 
r(a) 

where f(a) is a parameterization of the fixed worst trajectory, is a CLF for the system defined by (A, B). 
Hence, we have provided a 5-parametcr SSF in the (CC.2.2) case. The five parameters are: JC, the ratio 
Pb/ 'pa, and the three parameters involved in the change of coordinates to get the normal forms (15), (16). 

Remark 14 Notice that, in the cases (CC.l) and (CC.2.1) (cf. Section 3 and Theorem 3 in Appendix A.l), 
one can choose as SSF the set of quadratic polynomials, which actually is parameterized by two parameters. 

Remark 15 Let us come back to the general system (3), subject to HO. Notice that the question of finding 
the smallest m such that there exists a m-parameters finitc-SSF, for a certain class C of systems, has no real 
meaning if one does not require suitable conditions on the map \& in Definition 3. Indeed, it is always possible 
to build a countable SSF for the class of systems of type (3) in R" subject to HO. 



A Stability Conditions for bidimensional Systems 
A.l Statement of the Stability Conditions 

Theorem 3 Let A and B be two real matrices such that HI, H2, H3 and H4, given in Section 3, hold and 
define pa, Pb, JC,D as in Definition 4- We have the following stability conditions: 

Case (CC) If A and B have both complex eigenvalues then: 

Case (CC.l) ifV < then (V) is true; 
Case (CC.2) ifV>0 then: 

Case (CC.2.1) if JC < -1 then (V) is false; 

Case (CC.2.2) if JC > 1 then (V) is true if and only if it holds the following condition: 



Pcc : = ex P 



-pA arctan ( + pb j _ 



V 



,'pA-pBtC\ 7T 

Pb arctan [ — = — I - -{Pa + Pb) 



(PAPB + JC) + VV < l 

(paPb +fC)- ^/T> 
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Case (CC.3) IfV = then (V) holds true or false whether JC > 1 or K. < -1. 

Case (RC) If A and B have one of them complex and the other real eigenvalues, define \ '■= Pa^C — Pb, where 
PA and pb are chosen in such a way pa € «K, Pb G R- Then: 

Case (RC.l) ifV>0 then (V) is true; 

Case (RC.2) ifV<0 then x ^ and we /icwe: 

Case (RC.2.1) i/x > £/ien fPJ is false. Moreover in this case JC/i < 0; 

Case (RC.2. 2) if x < 0, then: 

Case (RC2.2.A) ifJC/i < i/ten (V) is true; 

Case (RC2.2.B) if JC/i > then (V) is true iff it holds the following condition: 



_l_ — ^(pA/i— 1) 

Pijc := (^r) e_PBi (Vl - ^ 2 ™~ sini - (cost - - sini)) < 1 (12) 



where: m — 



t = arccos 



{-PA/i ~ l)JC/i 

-PaH + PBlC/i 



Case (RC.3) IfV = then (V) holds true whether \ < or x > 0. 

Case (RR) J/ A and B have both real eigenvalues then: 

Case (RR.l) ifV<0 then (V) is true. Moreover we have \JC\ > 1; 
Case (RR.2) ifV>0 then JC ^ —PaPb (notice that —paPb > and : 

Case (RR.2.1) if JC > —paPb then (P) is false 

Case (RR.2. 2) if K < -paPb then: 

Case (RR.2.2.A) if JC > -1 then (P) is true; 

Case (RR.2.2.B) if JC < —1 then (P) is true iff the following condition holds: 

P R R ~ -} Sym {pA,PB,K,)f aSVm (pA,PB,JC) X (13) 

r ym (PB,PA,ic)<i, 

where: 



1 + pa/i + Pb/i + K, + \/V 
I pb/i - JCpA/i - \l> \ 



r™{PA, PB ,K) - 

f avm (pA,PB,JC) : 

\Pb/i - rCpA/i + W ) 

Case (RR.3) IfV = then (V) holds true or false whether JC < —paPb or JC > —paPb- 

Finally, if (V) does not hold true, then in case CC.2.2 with p cc = 1, case (RC.2.2.B), with P RC = 1, case 
(RR.2. 2. B), with p RR = \, case (CC.3) with JC < -1, case (RC.3) with x > and case (RR.3) with 
JC > —PaPb, the origin is just stable. In the other cases, the system is unstable. 

Remark 16 Formula (13) is a corrected version of Formula (6), p. 93, of [7] and it is proved in Appendix A. 3. 
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A. 2 Normal Forms of 2 x 2 Matrices 



Proposition 3 Let A, B be two 2x2 real matrices satisfying conditions HI, H2, H3 and H4 given in Section 
3. In the case in which one of the two matrices has real and the other non-real eigenvalues (i.e. the (RC) 
case ), assume that A is the one having real eigenvalues. Then there exists a 3-parameter change of coordinates 
and two constant a a, as > such that the matrices A/ a a and B/c*b (still denoted below by A and B) are in 
the following normal forms: 

Case in which A and B have both non-real eigenvalues ((CC) case): 

A -(-? B =(~? -»)• (14) 

where pa, Pb > 0, \E\ > 1. In this case, K, — h(E + -g). Moreover, the eigenvalues of A and B are 
respectively ~Pa i i and —ps ± i- 



Case in which A has real and B non-real eigenvalues ((RC) case): 

i/i + 1 
-p A /i - 1 



R _ / -p B-K,/ i -VI - /C 2 \ , , 

vT^ -PB + K/iJ' (16) 

where pb > 0, Pa/* > 1, /C 6 iR. In iftis case, i/ie eigenvalues of A and B are respectively —pA/i ± 1 
and — /9_b ± i ■ 

Case in which A and £> have both real eigenvalues ((RR) case): 

A - (~ PA o +1 <"> 

B - (\T/ (18) 

(19) 

where pA/i, Pb /i > 1 arl ^ /C G M \ {±1}- In this case, the eigenvalues of A and B are respectively 
—Pa/HI- and —pB/iH- 



A. 3 Proof of Formula (12) 

In this paragraph, we prove Formula (12), i.e. in the (RC.2.2.B) case, we determine an inequality defining 
the set of parameters pa, Pb, IC such that the property (V) , stated in Section 3, holds. 

Thanks to Proposition 3 (see also [7], Appendix B, p. 110), we can find a coordinate transformation such 
that (up to a rescaling of the matrices) A and B are given by equations (15), (16). In the case (RC.2.2.B), 
we have V := K? + 2p A pB^ - (1 + p\ + Pb) < °> X : = PA& - Pb < 0, JC/i > 0. Moreover, the set Q _1 (0) is 
the union of two lines passing from the origin and, at each point of Q _1 (0), the two vector fields point in the 
same direction. One easily checks that the slope of the two lines defining Q _1 (0) is: 

± -x±\/^P 
m = , 

c-pa/» - lyr^K* 

Notice that, in our case we have rri^ < and m + < m~ . 

In this case, the worst trajectories are concatenations of arcs of integral curves of the vector fields Ax, Bx 
and rotate counterclockwise around the origin. More precisely, they are integral curves of Ax from the line 
x-i = m + xi to the line X2 = m~x\, and integral curves of Bx otherwise. 
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Therefore, starting from the point ( ^ + ^ (with the field Ax), we follow the worst trajectory until it 

touches again the line X2 = m + x\. Property (V) is then satisfied if and only if p RC < 1, where p RC is the 
absolute value of the first coordinate of the final point. 

One can easily compute that the first switching time is t\ = — log . which is positive since — — > 1. 

2m' m 

Moreover, the integral curve of Bx starting from the point ( _ ^ is: 

-PBtf —\T\~~-~K? rn~ sin t + (cost — £ sint) 
V y/l - K? sint + m~(cos£ + ^ sini) 

and, setting the ratio between the second coordinate and the first one equal to m + , one obtains that the second 
switching time is t 2 = arccos — - . Notice that t 2 is well defined if and only if T> < (condition 

which is satisfied in our case). Moreover, i 2 is positive and less than n. Finally, the inequality we was seeking 
for is: 

/m +\~^ [pA/i ~ 1) / , JC \ 

Prc = \~) e~ PBt2 [yi-K?m- sini 2 - (cosi 2 - — sini 2 )J < 1. 

B The Stability Properties of (3) Only Depend on the Convex Hull 
of A 

We provide here the proof of Proposition 1. First, let us show the following: 

Claim Consider the switched system (3), under HO, and let A' be a measurable subset of A such that 
the convex hull of A' contains A. Then the following two conditions are equivalent: 

i) the system is GUES (rcsp. uniformly stable), with A(.) measurable, taking values in A, 

ii) the system is GUES (resp. uniformly stable), with A(.) measurable, taking values in A'. 

Proof of the Claim. Let A (rcsp. A') be the set of measurable functions ^4(.) : [0,oo[— * A (resp. A{.) : 
[0, oo[— > A'). Since A' is contained in A then the implication i) => ii) is obvious. 

Let us prove the other implication (that is strictly related to the classical approximability theorems in control 
theory) . We start considering uniform stability. By contradiction, assume that we can find e > satisfying the 
following. There exists a sequence of points (xi) tending to zero and a sequence of controls Ai{.) G Ji such that 
the corresponding trajectory 7; starting at xi exits the interior of the ball of radius e for some time t\. Using 
classical approximability results (see for instance [2]), the trajectory 7; can be approximated in the L°°-norm 
on [0, ti] by a trajectory 7 ; ' corresponding to a switching function A[{.) G Ji' and starting at x\. Hence 7^ exits 
the interior of the ball of radius e/2 at time ti. We reached a contradiction. 

Now we want to prove that GUES holds in the case A(.) £ Ji' implies GUES holds in the case A(.) £ Ji. 
Since A is compact, we know (see Definition 1 and below) that attractivity and GUES are equivalent for 
the corresponding switched system. Therefore, proceeding by contradiction, we can assume that there is a 
trajectory 7(.) of the switched system corresponding to A(.) £ Ji not converging to zero. That means that 
there exist e > and a sequence t n of times tending to infinity such that |7(t n )| > e. As before, we can 
approximate 7(.) on the interval [0, t n ] with a trajectory 7„(.) corresponding to controls taking values in A', in 
such a way that |7„(t„)| > e/2 . But this is impossible since we have assumed GUES for the switched system 
with A(.) £ A'. 

Notice that one can provide an alternative argument for the GUES part of Proposition 1, by using CLFs. 

□ 

Then one immediately extend to semicones observing that the stability properties of the system (3), subject 
to the compactness hypothesis HO, depend only on the shape of the trajectories and not on the way in which 
they are parameterized. 
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